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We give an exact, analytic, and manifestly gauge invariant account of pair production in combined 
longitudinal and transverse electromagnetic fields, both depending arbitrarily on Hghtfront time. 
The instantaneous, nonperturbative probability of pair creation is given explicitly along with the 
spectra of the final particle yield. Our results are relevant to high-intensity QED experiments now 
being planned for future optical and x-ray free electron lasers. 



I. INTRODUCTION 

The use of laser light sources in examining the high 
intensity regime of QED continues to draw attention, 
prompted by the advent of a new generation of both x-ray 
and optical laser facilities such as the European XFEL 
and the Extreme Light Infrastructure project [2 [2]. 
Dominating the theoretical activity in this area is the 
pursuit of (Schwinger) pair production in strong back- 
ground fields 

In this paper we apply the Dirac-Heisenberg-Wigner 
(DHW) formalism [3H5] to the problem of pair produc- 
tion in background fields. Within this approach, which 
was developed in the early 1990's [71 [S] and has gained 
increased attention in recent years [9Hl2]. one studies 
phase space distribution functions instead of the usual 
S-matrix elements. Since the DHW formalism deals es- 
sentially with quasi-probabilities, interpretation can be 
challenging. Nevertheless, quite some progress can be 
made by considering simple, but nontrivial, cases. 

An advantage of this approach is that the key object, 
called the DHW function, which is essentially the field 
density in an appropriate state, can in principle be ob- 
tained by solving a single partial differential equation. 
This is true at least when the gauge field is external: 
going beyond this is notoriously hard within the DHW 
approach |13| . while perturbation theory is obviously 
straightforward in field theory. In order to ensure that 
the obtained solution is physical, one must use sensible 
initial data which corresponds to calculating the DHW 
function in a physical state. Alternatively, one could try 
to construct the DHW function directly. This is chal- 
lenging since it requires finding the exact solutions of the 
Dirac equation in the chosen background field and then 
quantising the theory. 

In this paper we continue the investigation of the DHW 
function started in jTD] using hghtfront methods (see [TH - 
[T6] for applications of related methods to QED in a va- 
riety of strong external fields). Our previous results fo- 
cussed on plane wave backgrounds, i.e. transverse, or- 
thogonal electric and magnetic fields of equal magnitude. 
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depending on Hghtfront time a;+ . While we made progress 
in understanding the effective mass of a particle in an 
arbitrary pulse, we were of course unable to study pair 
production since single plane waves cannot produce pairs. 
In this paper we retain the plane wave fields but add a 
longitudinal electric field. Both of our fields will depend 
arbitrarily on hghtfront time, allowing us to model mod- 
ern short-duration laser pulses [TH [TTHTn] • 

To obtain the DHW function, we follow the second 
approach described above: we will therefore present 
new solutions of the Dirac equation in a combination 
of longitudinal and transverse fields, quantise the the- 
ory and calculate the hghtfront DHW function directly. 
This function, as we will show, can be interpreted as a 
(quasi)probabilistic measure of electron/positron occu- 
pation numbers. This will give us a clear signal of pair 
creation as we will be able to see, in a gauge invariant 
manner, the filling of states as particles are produced. 
We also confirm the results given in pUl - B5] , namely that 
from the infinitely boosted hghtfront frame one sees only 
the created positrons (modulo the choice of field) since 
the electrons decouple from the theory after creation. 

We begin in Sect. [H] by reviewing the DHW approach 
and presenting some basic results. In Sect. |III| we give 
the required solution of the Dirac equation in our chosen 
background. The quantisation of the theory and con- 
struction of the DHW function is not too hard but the 
expressions involved can become quite lengthy, and are 
therefore relegated to the appendices. We give the exact 
DHW function in Sect. |IV| and analyse pair creation on 
the hghtfront, using explicit examples of both short and 
long pulses. We also reconstruct the final particle spec- 
trum in the lab frame. Conclusions are given in Sect. |V] 
Our hghtfront conventions are quite standard but we en- 
courage the reader not familiar with hghtfront methods 
to consult the appendices for details. 



II. THE LIGHTFRONT DHW FUNCTION 

The DHW function is essentially given by the Fourier 
transform of the fermion field density in a chosen state, 
usually the vacuum. Our background fields will depend 
on hghtfront time a;+ = (x^ + x^)/V2 and so it is natu- 
ral to work in hghtfront field theory. On the hghtfront, 
the fermion fields ip split into dynamical fields "0+, and 
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constrained fields V'- (this is reviewed below, details are 
not needed here). It is therefore convenient, and simpler, 
to study the DHW function defined by the density of the 
dynamical fields rather than the full Dirac spinor. 

The equal lightfront time DHW function begins with 
the dynamical fermion density U in, say, the vacuum | ) . 
Noting that we use a sans-serif font to denote the spatial 
lightfront variables and momenta, i.e. x = {x^,x^} and 
p = this density is 

Uo.p = {0\[^^^ix+,X2),^Pl^ix\xi)]\0) , (1) 

where a and /? are spin indices. Setting X2 = x-|-y/2 and 
Xi = X — y/2, the DHW function is defined by Fourier 
transforming with respect to the relative co-ordinate y: 



j3 ^ip.y+ieJdz.A{z)jj 



(2) 



where gauge invariance of is ensured by the Wil- 
son line in the exponent. The line integral is taken over 
the straight path from xi to X2 which corresponds to 
minimally coupling the free DHW function by replac- 
ing 9 — >■ -D [21]. The factor of \/2 is a normalisation. 
Note that, being gauge invariant, the DHW function can 
be calculated in any gauge. Some intuition for what the 
DHW function represents can be built up by calculating 
it in a variety of simple states, to which we now turn. 



A. Free theory 

All our DHW functions will be proportional to the 
lightfront projector A+ = , so we write 



(3) 



The DHW function ^ for free fermions is easily found by 
writing down the mode expansion, calculating the expec- 
tation value and performing the Wigner transformation 
(without Wilson line). One finds 



W{x+;x, p) = Sign(p_) , 



(4) 



which is spatially homogeneous and displays only a sim- 
ple dependence on the lightfront momentum p_ . This 
behaviour is due to the existence of both positrons and 
electrons. To see why, replace |0) in ([T]) with |full), in 
which every positron and electron state is occupied. The 
DHW function becomes 
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FIG. 1. The free VF as a function of p_ , calculated in the ordi- 
nary, empty vacuum (black), in the state filled with positrons 
(red, dashed), the state filled with electrons (red, dotted) and 
the completely filled state (blue). The region p_ < to the 
left of the plot is controlled by positrons, p_ > by electrons. 
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FIG. 2. The DHW function in the mixed state, Q. The 
deviation from the vacuum DHW function is given by twice 
the probability for being in the state j full). 



is a matter of convention (it does not refer to negative 
energy) which follows from the choice of exponent in the 
transform ([2]), since the mode expansion for -0 looks like 



Following this, the DHW variable p may be associated 
with the momentum of an electron, or minus the mo- 
mentum of a positron. Consider also a mixed state. 



mixed ) = \/l - PI ) + VpI full ; 



(8) 



which has probability P of being full and 1 — P of being 
empty. The DHW function is easily found to be 



W{x^;x, p) = Sign(— p_ 



filled vacuum 



(5) 



Similarly, filling all the electron or positron states, one 
obtains instead 



VF(a;+;x, p) = -1 

= +1 



electrons filled , 
positrons filled 



(6) 



These results are shown in Fig. [T] The region p_ > 
is controlled by electrons and p_ < by positrons: this 



W = Sign(p_) + 2PSign(-p_) 



(9) 



so that, for example, the p_ < portion of the DHW 
function is raised by a height 2P relative to that in the 
vacuum, see Fig.|2] The DHW function therefore gives us 
a (quasi)probabilistic measure of the occupation numbers 
of electrons and positrons. The above points are worth 
keeping in mind for later, as they will aid our interpreta- 
tion of the DHW function in the interacting theory. 
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B. Plane wave backgrounds 

The DHW function for both scalars and spinors in an 
arbitrary plane wave background was calculated in |10j . 
That paper also considered the covariant DHW func- 
tion in which the lightfront times are also separated, 
i.e. one works with the density [tp+{x2),'4'\{xi)\ with 
x^'^a = 2;'' ± 2/''/2. The DHW function is then defined 
by an integral over d'^y, and was found to be 



Sign(p_) / dy+ exp 
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PI 



2p- 



(10) 

where Af ^ is clearly an effective mass (see [lOl |25] for de- 
tails) which extends the intensity-dependent mass shift 
from purely periodic plane waves |26H28| . to arbitrary 
plane wave backgrounds. Integrating out one recov- 
ers, by definition, our W: this is precisely the same as 
in the free theory Q. The natural interpretation of this 
result is the well known statement that plane waves do 
not create pairs. 

In the following sections we will examine what happens 
to the DHW function when a pair-creating longitudinal 
electric field is added to the plane waves. We will con- 
tinue to focus on the simpler W rather than the covariant 
DHW function as this is the more common approach in 
the literature, and because, following the above, any de- 
viation in W from the free vacuum result Q must be 
due (at least in part) to the longitudinal field. 



III. LONGITUDINAL AND TRANSVERSE 
FIELDS 

We consider a laser pulse moving up the x^-axis. This 
defines our 'longitudinal' direction, while x^,x'^ are the 
transverse directions. A variety of models for the laser 
field can be found in the literature: the case of a con- 
stant, longitudinal electric field is of course covered by 
Schwinger's classic result 29J. The models for which 
most analytic progress can be made (in terms of calculat- 
ing scattering amplitudes) are plane waves depending on 
X'*' [2SH1H]. The combination of a constant longitudinal 
electric field and periodic plane waves was described in 
P5] . The case of a purely longitudinal electric field E[x*) 
depending arbitrarily on x^ was covered by [501 HI] (and 
includes Schwinger's result as a particular case). The 
lightfront methods used in those papers can be extended 
to cover the case of longitudinal E{x*) with a parallel 
B-field also depending on x^ ^30j. For fields depending 
on both a;+ and x^ see [3T1[32]. Purely longitudinal elec- 
tric fields E{x^) depending arbitrarily on (instant) time 
x^ are widely used in the literature. Such fields model 
the focus of counter propagating laser pulses in which the 
magnetic field components cancel. They have been used 
to investigate pair production for oscillating fields [33l - 
135] , pulsed fields [SB] [SZ] and pulsed fields with sub-cycle 
structure 



Here we further extend the above results, covering the 
case of longitudinal electric and transverse electromag- 
netic plane wave fields, both depending arbitrarily on . 
We work in 'anti-lightcone' gauge A~ = = (the 
usual lightcone gauge is A_ — 0). The remaining com- 
ponents of the potential are given by. 



A_=-jdyE\y), A^ = V2 J dy E^ {y) , {U) 



where we assume for simplicity that the fields turn on at 
x+ = 0. This can, and will, be relaxed below. 



A. Solutions of the Dirac equation 

Defining the projectors A* = ^7^7* the fermion field 
decomposes into ■0 = V'+ + V'- with ip± = A^ip. The 
Dirac equation then separates into 

id+i;+ = \{i-i^D^ + m)r^- , (12) 
iD_'4}_ = \{i"f^D^ + m)7+V+ ■ (13) 

We immediately take the Fourier transform of the trans- 
verse coordinates, id^ — > k^, which replaces 



eA^{x^) EE 7r^(2;+) 



(14) 



One solves the Dirac equation by first observing that ■0- 
is a constrained field, since it can be expressed in terms 
of 1/'+ using ( 12 ): 



_ 7^7r_L -f 771 . 

tp- = 5 iTo+ip+ 



where the mode frequency is defined by 

Lo\x^) = Trl{x+) + 

= kl+ ir? + e^A\{x^) - 2ek^A^{x^) 



(15) 



(16) 



in which we recognise the Volkov exponent |41| . Substi- 
tuting (151 into (131, and noting that D_ and 9+ do not 



commute, one obtains a simple equation for 

(17) 

we see that 



If we try to Fourier transform i9_ 
solving (17 1 requires inverting 



k_ 



k_ - e^_(a;+) 



(18) 



which can clearly vanish, possibly multiple times, for a 
given k_. This is the zero- mode problem of lightfront 
field theory [32], but made time- dependent by the exter- 
nal field. The physics of the zero-mode in the current 
context is as follows. An electron with momentum k_ at 
time x^ = Q acquires (as follows from solving the Lorentz 
equation) a momentum k'_ = k_ — eA_{x'^) at later times. 
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from (IT5I). The function G is 



FIG. 3. The domain of our solution (19l, > and x > 



~L. The finite-duration baclcground fields depend on x G 
0...x^, as is also illustrated, along with the behaviour of 
electrons and positrons created within the field. 



We have e < 0, so if we imagine that E^^ is positive, then 
k'_ will vanish at some later lightfront time, which means 
the electron moves parallel to the axis (reaches the 
speed of light) at this instant: it therefore vanishes from 
the theory since it cannot be seen at any subsequent light- 
front time [22 [23] . Note that a positron's momentum, 
on the other hand, only increases in the above circum- 
stances. Hence the positrons remain in the theory. This 
will be useful for later. 

We are now ready to give the solution of the Dirac 
equation. We follow the method of [201 l2l]- The idea 
is to turn the fields on at a;+ = (for convenience, this 
can be relaxed, see below), and solve the Dirac equation 
in the semi-infinite region > 0, x~ > —L for some 
positive L, in terms of initial data. This gives, as we 
will see, a prescription for handling the singularity at 
k_ — eA_ = 0. In the end the limit L — > 00 is taken. 
The domain of our solution, together with an illustration 
of our fields and the motion of particles within them, is 
shown in Fig. [3] 



It may be checked directly that the solution to ( 12 ) 
pl| in x+ > 



> -L is 



V'+(x+,x-) = / dy- ^P40,y-)D4y)G{x\0;x-,y-) 



dy+ 



{y\-L)G{x\y^;x-,-L) , (19) 



with given by (15). We consider the various terms. 
First, V'+'s dependence on the boundary data is explicit: 
the solution depends on on the characteristic a;+ = 
and d+ijj+ ^ "0- on the characteristic x^ ~ —L, since 

id+^lj+ = lif + m.)j-^_, (20) 



G{x+,y+;x ,y ) 



J(V -X )(k_+i/L) 



27r k_-eA_{y+) 
and £ is defined by 



(21) 



£k-{x'^,y'^) = exp 



i/L 



ds u}^{s) 
- eA_(s) + i/L 



(22) 



It is worth considering this function in a little detail, as 
it exhibits the essential difference between the transverse 
and longitudinal fields. The transverse plane wave fields 
enter just as in the Volkov solution, in the numerator of 
the exponent [41^. These terms may therefore be recov- 
ered by resumming all orders of perturbation theory in 
the plane wave coupling eAj^. The longitudinal field, on 
the other hand, appears in the denominator and exhibits 
a singularity on the real line, regulated by the factors of 
i/L: when L — 00 this leads to an essential singularity 
in the coupling, as in Schwinger's result. 

An advantage of the approach we adopt is that the 
differences between the types of field, and the important 
structures, are laid bare. Nothing is hidden inside the 
behaviour of special functions, as is frequently the case 
in the instant-form approach: the equal x'^ (instant time) 
DHW function is expressed in terms of parabolic cylin- 
der functions for i?(x°) = E, constant, and in terms of 
hypergeometric functions for E{x^) = Esech^{ujx^) [5]. 
However, a disadvantage of our approach is that expres- 
sions quickly become lengthy. For this reason, the quan- 
tisation of ( [Tq] ) and the calculation of the DHW function 
are left to the appendix. Related calculations are explic- 
itly performed in .21j . which the reader may consult for 
further examples. The final result for the DHW function 
in the limit L — c» is, however, extremely compact, and 
we turn to it now. 



IV. PAIR PRODUCTION 

If one considers only a longitudinal electric field de- 
pending on one finds that not only the vacuum per- 
sistence amplitude but also the pair production rate may 
be calculated instantaneously as a function of a;+. The 
derivation of this latter result requires a careful interpre- 
tation of the Heisenberg operators in order to identify 
the pair production probability [201 [H] . We can provide 
a (positive) check of that interpretation using the DHW 
function. Moreover, our approach is manifestly gauge 
invariant. 

It was found in |28| that the addition of a plane wave 
to a constant electric field (Schwinger's case) does not 
change the vacuum persistence amplitude. We will see 
for our fields that this remains true: the plane wave has 
no impact on the creation of particles, nor the properties 
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they are created with. Rather neatly, though, our results 
make clear the post-creation effect of the plane waves on 
the particles. 



A. The DHW function 

Our solution of the Dirac equation (and its quantisa- 
tion, as described in the appendix) is valid for arbitrarily 
longitudinal and plane wave fields. We present here re- 
sults for the case in which i?" is assumed to be positive, 
so that eA_ is positive and increasing, as this is when the 
singularities in (22 1 have the simplest structure in phase 
space. Our longitudinal fields therefore model subcycle 
pulses, which are of considerable current interest |43j . 
The plane wave fields remain arbitrary. 



We can now give the DHW function W. At 



0, 



W is that of the free theory, see (|4| . Once the longitudi- 



nal fields turn off at 



the theory becomes stable 



against pair production and we find that the DHW func- 
tion again becomes constant in lightfront time, matching 
its final value in the pulse. For the duration of the pulse, 
the DHW function is 



i.e. < x+ < 



W(a;+;x,p) ^ Sigii{p_) + 2F 9{-p_)e{eA_{x+) + p_) . 

(23) 

The first term is the DHW function of the empty vacuum. 
The second term contains the effects of the background 
fields, and has a form similar to that in ([9|, though it is 
restricted to p_ < 0. The description and investigation 
of this term, in particular P, will occupy the remainder 
of the paper. 

It is important before embarking on this to give the 
interpretation of the DHW variable p. In the free the- 
ory, see Sect. |ll] p (— p) is the kinetic momentum of an 
electron (positron). This also holds once the background 
fields turn off and the theory again becomes free (the fi- 
nal particle spectrum is of course what we would be inter- 
ested in experimentally). The DHW function smoothly 
connects the initial and final distributions in a gauge in- 
variant manner. Furthermore, ^ shows that (canoni- 
cal!) momentum dependence on, say, k in the density 
will be set equal to p + eA in the Wigner function, so 
that p is naturally interpreted as a kinematic momen- 
tum. From here on we therefore associate p_ < with a 
positron momentum 7r_ = — p_ > 0, as in the free theory: 



from (23), this is clearly the region of interest. 



B. From dynamics to probabilities 



initial time x^. From the Lorentz equation, it will at a 
later time a;+ have momentum 



= eA_{x+) - eA_{x^) 



(24) 



using (11). Since eA_ is positive and increasing for the 
duration of the pulse, it has a unique inverse Xp such 
that eA_{Xp) — p and A^eA_(2:) = It follows that, on 
observation of a positron with momentum tt_ at time x^ , 



the 'initial time' could be reconstructed from (24) 



eA_{x^) = eA_{x+) - 7r_ 



(25) 



(We suppress the dependence of x^ on x"*" for compact- 
ness.) If the positron also has zero transverse momentum 
at the initial time, its later transverse momentum is 



7r_L(7r_) = eyl^(x+) - eA^(a;+(7r_)) 



(26) 



(Again suppressing dependence of tt^ on a;+.) It is 



clear from the integral expressions (11) that the results 



(24)-(26) simply describe the energy transferred to the 



positron from the background fields over the elapsed 
time. This is illustrated in Fig. |4] 

With these definitions we can give a very simple ex- 
pression for P: 



exp 



TTTO^ +Tt[p^ + 7r^(— p_)]' 

|e|Ell(x+(-p_)) 



(27) 



where x^ and ttj^ are defined in ( 25 ) and ( 26 ) . We recog- 



nise a similar structure as found in Schwinger's results, 
but for more general fields, and also depending instan- 
taneously on lightfront time. Comparing (23) and (|9|, 
P is naturally interpreted as the probability that the 
positron states with momentum 7r_ = —p_ have been 
filled by time since particles are being created by the 
background fields. Using the dynamics discussed above. 








{0,0} 



{7r_L , 7r_} 



In order to give the most compact and intutive expres- 
sion for P it is useful to recall some results on the motion 
of particles in our background fields. The Lorentz equa- 
tion for a positron with kinematic momentum tt^ and 
charge — e > is dvr^ = —eF^^dx'^. Suppose then, that 
a positron is created with momentum 7r_ = at some 



FIG. 4. A particle is created at time 3;+, with probability 
determined by the electric field strength at that time. The 
particle has zero longitudinal momentum, and transverse mo- 
mentum normally distributed around zero. It is observed at 
a later time , after which it has acquired longitudinal and 
transverse momenta 7r_ and 7r±. 
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a more precise statement is the following: P gives the 
probability of observing positrons with momenta 7r„ and 
TT^ at time such that these positrons were created 
at with 7r_ = and transverse momentum normally 
distributed about tt^ = 0. Note that because the argu- 
ment of the pair creating field in P is not x+ but , the 
probability of observing a positron is dependent on the 
electric field strength at the moment of creation x^ , and 
not on the 'observation' time see also Fig. |4] This is 
a neat and physically sensible result. 

We now explain this interpretation in more detail and 
reinforce it with a series of examples. In order to keep the 
presentation as clear as possible, we begin by dropping 
down to 1+1 dimensions, turning off all transverse depen- 
dence. This will be reinstated below. The first obvious 
question, and obvious difference between (23) and (|9]), is 
why do we not see electron states being filled? Why is 
there no change to the DHW function for p_ > 0? 

A related result was found and explained in [jQl HI] , 
which we now describe in our language. Our DHW func- 
tion has the form of a disturbance propagating out from 
p_ = as time evolves. This is also where the singu- 
larities in the Dirac equation live in the L ^ oo limit, 
and we have already seen that P describes particles with 
zero initial longitudinal momentum. Put together, this 
means that pairs are created travelling at the speed of 
light. The distinction is that the electrons, being acceler- 
ated down the x'^-axis by the positive field, travel parallel 
to the x~ direction and so, from the perspective of the 
infinitely boosted lightfront frame, immediately vanish. 
The positrons, on the other hand, are accelerated u p th e 



-axis and therefore acquire positive 7r_, see also (24), 



and remain visible in the lightfront frame. We therefore 
confirm the result of [20, 21 : from the lightfront perspec- 
tive, one only sees the positrons. 



The final piece of (23) to consider is the second step 
function. This states simply that the argument of E\ 
that is ^eA(2;+)-7r_ ; must be positive; in other words, 
sufficient time must have elapsed for a particle with zero 
longitudinal momentum to have acquired 7r_ by time x^ . 
Together, the two theta functions therefore imply, using 
p4|-(l26|, that 



< x+ < a;^ 



(28) 



which is a simple statement of causality: observed pairs 
must have been created at earlier times, but after the 
longitudinal field turns on. We conclude that the DHW 
function shows us pair production from the vacuum, in 
real lightfront time, with P the probability of pair cre- 
ation. The positrons appearing at time x^ with a fi- 
nite range of momenta 7r_ are subject to the (natural) 
constraint that sufficient time must have elapsed for the 
positron to have absorbed this momenta from the fields, 
starting from 7r_ = 0. We now move on to explicit ex- 
amples, staying in 1 4- 1 dimensions for the moment. 




FIG. 5. The DHW function in the subcycle pulse 1 29 1, plot 



ted as a function of p_ , for zero transverse momentum. The 
maximum allowed momentum, see ( |31| ) is also shown in the 
(p_ , i^) plane. 



C. Example: finite pulse duration 



We begin with the electric field 

E\x+) = Eosm{ujx^) 



(29) 



for < ujx'^ < TT and zero otherwise, modelling a half- 
cycle of the laser. The DHW function W is plotted in 
Fig. [5] When plotting, we use rescaled variables 



\e\Eo 

m 2 



P- 



p_uj 



(30) 



which measures the electric field strength in units of the 
Schwinger field and p_ in units of (as we are about to 
see) half its maximum value. We have chosen = t: 
for our plot, which means our electric field strength is 
roughly three times higher than the Schwinger limit: this 
compensates fully for the damping factor in the exponent 
of (27) and allows us to clearly see the behaviour of the 



DHW function. We consider other field strengths below. 

In Fig.jSj we see that W matches the free theory result 
([4]) at = 0. As time evolves, W becomes both p__ and 
x+ dependent, with the deviation from vacuum spreading 
out from p_ = at x~^ = 0. At time x^ the function is 



explicitly limited in extent by the theta- functions in ( 23 ) , 
which give 



eA_{x+) +p_ >0 



> X. 



(31) 



The DHW function eventually stabilises as the field 
switches off at x'^ ~ ^'f^ upon which it is straightforward 
to extract properties of the final positron distribution. 
The final range of possible positron momenta is dictated 
by, following the above, —p^ > and —p_ < eA_{x'j:). 
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FIG. 6. The DHW function in the sech^ pulse |37| with 
peak amphtude eo = 1, see (30 1, plotted as a function of p_, 
zero transverse mom entum. As before, the range of allowed 
momentum, see (31 1 is shown in the (p_,£+) plane: there is 



a smoother falloff than in the previous example. 



D. Example: adiabatic switching 

In the light of recent literature results, which we dis- 
cuss below, it is worthwhile pointing out that there 
is nothing to stop us turning our fields on arbitrarily 
smoothly starting from arbitrary initial times, without 



affecting the essential properties of our solutions (19) or 



our DHW function ( 23 ): our choice of switching the fields 
on at x+ = was for convenience. We therefore consider 
a field which, while qualitatively similar to our previous 
example, falls off quickly but smoothly at ±oo, namely 



(37) 



where Eq gives the peak intensity and 1/uj the effective 
duration of the pulse. The corresponding gauge potential 
is given by the integral from a;+ = —cxd of this function, 
and is therefore 



eA_{x+) = ^^(l+tanh(wx+)) , 



(38) 



In our current example, lox^ = tt, so the final range of 
positron momenta is 



< 7r_ < 2^^^^ = 2mao , 



(32) 



where, in the final equality, we have introduced the peak 
field intensity oq |i5] . 



\e\E,, 



ao 



(33) 



The momentum distribution is peaked around that value 
of 7r_ such that the electric field is maximal at the instant 
of creation. Let Xq be the time at which i?" (xq ) is max- 
imal, then the most probable kinetic momentum {tt_) is 
the solution of the equation 

For our current example ujxq — tt/2 and we find 



(7r_) = = mao , 



(35) 



which is also clearly seen from Fig. [5] Since the momenta 
are on-shell in the free theory, we have (tt^) = TO^/2(7r_) 
and so we can easily convert these expressions back to 
cartesian co-ordinates to find the likely energy and z- 
component of the momentum, which are 



ZOq ZOq 



map . 



(36) 



(Note that the probability P for producing very high en- 
ergy particles 7r_ ~ in the range (32) is almost zero.) 



where the constant term follows from the definition (11 1, 
with the initial time translated to a;+ = — oo. The result- 
ing DHW function is plotted in Fig. [6] for a peak field 
strength equal to the Schwinger field. The DHW func- 
tion is shown over the temporal range —1/2 < ujx^ < 2. 
Before this, the deviation from vacuum is minimal and 
afterwards the DHW function becomes effectively sta- 
ble. The final range of positron momenta which may be 
created by the field remains finite and, with our choice 
of variables, is expressed just as before, < 7r_ < 2mao. 
The probability P is peaked, in the limit — >■ cx), around 
7r_ = mflQ. Note that even at the Schwinger limit, the 
probability of pair production remains small. 



E. Example: ever-increasing field 



We note that the factor of 2 in ( 23 ) can be understood 
from a second perspective. Suppose we consider an elec- 
tric field which increases without bound, for example 



E^^{x+) = EaLux+ . 



(39) 



As time increases, the electric field becomes overcritical 
and we expect the probability of creating particles with 
any given momentum (in the allowed range, which also 
expands in time) to approach unity. This means that the 
state should become filled with positrons, and we expect 
to recover, from ([6]), W — +1 in the region p_ < 0. This 
is precisely what the factor of 2 ensures: if P — > 1, the 
DHW function approaches W — 1 for p_ < 0. This is 
shown explicitly in Fig. [Tj the DHW function transfers 
from —1 to +1 over the whole negative p_ range as time 
evolves. 



FIG. 7. DHW function in the field ( 39 1, with eo = 0.5, which 



increases without bound. As time evolves, all positron states 
are filled with unit probability. 



F. Comparison with the instant-form approach 

Our results share some similarities with investigations 
of particle creation in x'^-dependent electric fields, within 
the usual (instant-form) DHW formalism. In that ap- 
proach, a certain combination of instant-form DHW 
spinor components can be interpreted as a particle num- 
ber density f44], and its behaviour is as follows. As the 
electric field grows with time energy is transferred to 
the Dirac field such that a peak around = develops. 
At intermediate (non-asymptotic) times this peak is in- 
terpreted as being composed of virtual electron-positron 
pairs. As time evolves, only a part of these virtual par- 
ticles become real particles which are then accelerated 
by the electric field and spread out from = 0. At 
asymptotic times a;*^ — > cxd, the particle number density 
of real particles stabilises whereas the virtual electron- 
positron peak around = disappears again. An ex- 
ample is shown in Fig. [8] for the asymptotically switched 
field E{x^) — _Eo sech^(a;a;°). For more details see [44] . 

One difference between these results and our own is 
that the lightfront DHW function does not exhibit the 
intermediate virtual particle peak or oscillatory struc- 
ture seen in the instant form. The reason for this seems 
to be that pair production on the lightfront is an in- 
stantaneous event, occurring at the instant when a given 
fermion mode can produce a particle of zero longitudinal 
momentum, see (18 1. This is confirmed by our expres- 
sion for the pair creation probability P: it is expressed 



entirely in terms of classical particle trajectories. 

This is quite intriguing, as it may be related to the 
'triviality' of the lightfront vacuum. The instant-form 
vacuum is filled with virtual pairs which can be pulled 
onto the mass-shell by the external field. Recall that 
the Schwinger field strength can be obtained by equating 
the electron rest mass with the work done by the elec- 



-1 20 



FIG. 8. Particle number density M for the x^-dependent 
sech'^ pulse within the instant-form DHW formalism, plotted 
as a function of ps and zero transverse momentum. Peak 
amplitude eo = 1- 



trie field over the lifetime of a virtual pair: in this sense, 
there is a time scale involved in Schwinger pair produc- 
tion. The lightfront vacuum, on the other hand, is often 
referred to as 'trivial', which is the statement that it is 
completely empty of particles, both real and virtual |46) . 
Moreover, it is stable. In this picture, then, pairs are cre- 
ated from the energy pumped into the system, not from 
virtual particles being pulled on-shell, and the Schwinger 
'time-scale' is absent. This is an investigation for another 
time, though. We now return to properties of the DHW 
function. 



G. Nonperturbative dependence 

Recent results on pair production in ^''-dependent elec- 
tric fields (in the usual instant time DHW formalism), 
find a purely perturbative dependence on a particular 
electric field which switches on adiabatically in the infi- 
nite past [17]. Moreover, it is stated that the essential 
singularity of Schwinger's results must therefore be due 
to the unphysical nature of a constant, ever-present elec- 
tric field. 

Let us reconsider our results in this light. We have seen 
explicitly that for both sharply and smoothly switched 
fields, the pair production probability is basically de- 
scribed by a factor exp(— m7r/|e|£'"), which retains the 
essential singularity in the coupling from Schwinger's re- 
sult. It is therefore clear that the nonperturbative nature 
of pair creation is not due to some unphysical assump- 
tion about when, or how smoothly, the fields turn on or 
off. The dependence of our results on e is of course a 
little more complex than that, since i?" is evaluated at 



4(-p_), see (27) 



It is useful to examine the form of the (final) probabil- 
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ity P when the fields turn off. We work with the physical 
momentum 7r_ here, and turn off the plane wave fields 
for simplicity, setting also tt^ = 0. We begin with the 
field (29), for which the final probability is 



loe 



7r_ (2mao — 7r_ ) 



(40) 



Note that the denominator is positive because of the fi- 
nite allowed 7r_ range. We can examine this probability 
for, for example, small momenta (which corresponds to 
extremely energetic particles in the lab frame) by expand- 
ing in 7r_ : 



1 



4|e|^;o 



(41) 



We clearly see the l/|e|i?o dependence. What happens 
for the adiabatically switched sech^ pulse? For this field, 
the final probability is 



■K\e\EQW? 

lim — log P = — — 

i:+->co 2|e|a;£'o7r_ — w 



(42) 



Again, the denominator is positive, and we can make the 
same small 7r_ expansion as above, finding 



lim — log 

r+ — ^Dn 



2w7r_ 



1 



2e£;, 







(43) 



This displays a different dependence on the various pa- 
rameters. In particular, the dominant term is indepen- 
dent of the field strength E^. Does this correspond to a 
perturbative dependence on the field strength? The an- 
swer is no: not only does the second term in the expan- 
sion contain explicitly nonperturbative (Schwinger-like) 
terms, but the leading term of (431 actually contains a 
hidden nonperturbative dependence. To see this, note 
that the leading term, despite not being explicitly de- 
pendent on _Bo, does not survive the limit i?o — > (which 
would contradict the free theory). This is because the 
range of ■n_ is finite, being limited by eA_^ so that tak- 
ing i?o — >■ implies taking 7r_ — >■ 0, and this kills P as the 
fields turn off. 

What this result really shows is only that, and as is 
not surprising, the distribution of the produced parti- 
cles depends on the geometry of the field, for example 
whether the field turns off sharply or smoothly. We have 
not been able to identify a regime where the results may 
be expressed as a perturbation in Et^. We stress that this 
holds at least on the lightfront, for fields depending on 
x'^: there are differences between this and the instant- 
form approach, see Sect. |IVF[ above. 

It is, though, entirely possible for the effective action to 
exhibit both perturbative and Schwinger-like nonpertur- 
bative behaviour when the electric field depends on 
(and is even adiabatically switched), depending on the 
relative sizes of the parameters involved. This is shown 
in (48j . which discusses many deep connections between 



perturbative and nonperturbative physics. One is lead 
to conclude that the perturbative dependence found in 
|47) , while very interesting, is not inconsistent with the 
existence of nonperturbative behaviour. Combining this 
with our own results, we do not believe that any doubt 
is cast on the validity of Schwinger's result. 



H. Transverse dependence 

Finally, it is time to return to 3-1- 1 dimensions proper, 
and allow for plane wave fields. Consider the full ex- 



pression (27) for the probability. The plane wave contri- 



butions do not appear in the step functions, thus they 
do not affect the constraints dictating the momentum 
ranges. Nor do the plane waves enter into the argument 
of the pair creating field £'". 

It is clear that without the plane waves, the probability 
for production of pairs with nonzero transverse momenta 
is normally distributed around tt^ = 0. Turning on the 
plane waves, it may seem strange at first glance that the 
peak of this distribution is shifted to nonzero values: is 
the plane wave affecting the probability of pair produc- 
tion? The answer is no: recalling (24)-(26), the distribu- 



tion in (27) is obtained for positrons created with trans- 



verse momentum normally distributed around tt^ = 0, 
and which at the subsequent time a;+ must have acquired 
transverse momentum 



X ■ 

y ds \/2e£;^(s) 



(44) 



from the plane wave fields, using the Lorentz equation. 
The DHW function therefore takes into account both 
what happens at the instant of creation, but also what 
would subsequently be observed. 

To summarise, the plane waves do of course influence 
the particles after they are created, and so it is no sur- 
prise that they appear in the DHW function and the 
final particle distribution. The plane waves do have no 
influence, though, on whether particles are created or 
not. This is reaffirmed by integrating over momenta to 
obtain, for example, the total probability of pair produc- 
tion or the vacuum persistence amplitude: one finds that 
all dependence on the plane wave fields vanishes because 
7r^(— p_) can be absorbed into p± by a change of vari- 
able. Thus the plane wave fields, and in particular any 
effective mass they may generate, do not influence the 
probability of vacuum decay [501 1211 HI] • 



V. CONCLUSIONS 

We have investigated the phenomenon of non- 
perturbative pair creation in background electromagnetic 
fields, within the lightfront DHW formalism. We calcu- 
lated the DHW function by solving the Dirac equation 
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in a combination of longitudinal electric and transverse 
plane wave fields which both depend arbitrarily on light- 
front time. This extends the work of [201 IHl 121] to an 
even wider class of fields. 

As shown in |10j, the DHW function W is not altered 
by a single plane wave field, since a plane wave can not 
produce pairs. Switching on an additional pair-creating 
electric field, however, one observes a deviation from the 
vacuum result which signals pair creation. The pair cre- 
ation probability itself is exponentially suppressed by a 
factor exp(— TO7r/|e|i?") which retains the essential sin- 
gularity in the electromagnetic coupling e in Schwinger's 
result, but is valid for much more general fields, and 
in particular is independent of how the electric field is 
switched on and off. This may be contrasted with recent 
results in the instant-form approach |47j . 

Notably, we have seen that the value of the DHW func- 
tion W is altered only in its positron sector, whereas the 
electron content remains unchanged. This confirms pre- 
vious results that, from the lightfront perspective, only 
one of the particle species remains in the theory following 
creation ^Ull^ . 

All in all, the DHW function can be a powerful 
tool in analysing quantum physics in background fields 
(particularly pair production), and in a language which 
is essentially classical. We have given an elegant and 
physical interpretation in terms of the pair creation 
probability and the subsequent dynamics of the parti- 
cles. This makes it clear that we observe real particle 
production in the lightfront formalism. 
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Appendix A: Notation. 

Our lightcone directions are = (x'^ + )/V2. We 
prefer momenta to carry covariant indices and so p_ is 
a spatial momentum conjugate to a;", while p+ is the 
lightfront energy, conjugate to lightfront time x+. The 
metric has determinant —1 as in Cartesian coordinates 
and contra-/covariant indices are related by = v±- for 
arbitrary vectors v. We use a sans-serif font to denote 
the 'spatial' variables and momenta, i.e. x = {x~,x^} 
and p = {p-,p±}- Our Fourier conventions are 

f{p) = Jdx e^^^fix) , fix) = e-'^V(p) , (Al) 

and we Fourier transform before taking conjugates so 

Pip) = J dx e-^P^Pix) , pix)^J^ e^^-pip) . 

(A2) 



Appendix B: The quantum theory 



We now wish to quantise our solution ( 19 1. Quantisation is performed by imposing canonical commutation relations 
on the initial data. As in [Ml HI] one finds that the quantised spinor fields obey the desired commutation relations, 
returning briefly to full co-ordinate space. 



{^Ax),i^liy)]l^^y, = ^Kn\x^-y-)5ix--y-) , (Bl) 
provided that the initial data obeys 



Y (B2) 

{^_(^),^l(y)}L-=,-=_i = -^^-5\x--y-)5ix^~v^) , 

What we are really doing here is solving a Cauchy problem with initial data on two lightlike characteristics: this is not 
quite what one usually does on the lightfront (where the operators on the x^ characteristic are not needed explicitly) 
but is necessitated by the time-dependence introduced into the zero-mode problem. Note that without the operators 
on x^ — —L one does not recover known results such as Schwinger's vacuum persistence amplitude in the L ^ oo 
limit. Moreover, the approach used here has been verified by alternative methods [5T]. (The question of whether 
there is a method to recover the results of this approach in ordinary lightfront quantisation without the operators at 
x^ — has not, to our knowledge, been addressed.) 

From here one can construct the Hamiltonian and states of the theory as normal. Computationally, it is useful to 



work in what becomes Fourier space in the L ^ oo limit: considering the A:_ -integrations inside the functions G, (21 1, 
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we define the 'almost' Fourier transform -0+ by 



M^\x-) EE / — . (B3) 

ZTT 



Dependence on k± will not be written unless it is needed. Explicitly, our new field is 



x+ 







fc_ - eA_{y+) + i/L 



(The operators become canonically normalised creation (and annihilation) operators as L — >■ oo, in the interpreta- 
tion of [2Q1I21], up to a factor of 2^/**: hence the factor of V2 in our DHW function The state of interest for us is 
the vacuum, or rather the free vacuum, since what we wish to do is begin in this state at — 0, apply our external 
fields and see what happens as time evolves. Since we are working in the Heisenberg picture, in which the states are 
time- independent, the state we need is precisely the free vacuum state. This means that the density we need, ([!]), is 
calculated by first expressing in terms of the free initial data using (B4), and then evaluating the density of these 
free fields in the free vacuum state. This calculation is performed using the ordinary free-field mode expansion. For 
example, a free-field calculation easily gives us 

oo 

1 . , rdk_ 



which we will use below. 



Appendix C: Calculating the DHW function 



We now wish to calculate ([T]) and To do so we make an assumption on the longitudinal electric field, namely 
that it is positive, e.g. modelling a subcycle pulse. As a result, eA_ is a positive function increasing from at 
x'^ — 0. This approximation can be relaxed, but doing so means that the analytic structure in the functions £ 
becomes significantly more complex, since the kinematic momentum 'k_{x'^) = k_ — eA_(a;+) may have multiple zeros. 
The 'subcycle assumption' means that this function has at most one zero for a given k_. With this assumption, we 
can define the inverse function X{k_) = Xk_ by 

k_''eA_{Xk_) = . (CI) 

This function exists for the duration of the longitudinal field i?". What happens after the fields turn off is explained 
below. To proceed, we express the density U in terms of the modes '0+(x''', fc_). The first term of the commutator is 



[0\i;4x\x^)i>lix\xl)\0 



dL 
2^ 



fil-)e 



2tt 



riq- 



Q\^^{x\l_)^X{x\q_)\Q) , (C2) 



where we have included test functions under the momentum integrals in order to take the L oo limit more 
rigourously. Since '0+ itself contains two terms, the overlap in (C2| contains four. Only one of these, the product of 
the first term in (B4| with its conjugate, ultimately contributes, so we present only this calculation. Using (B4| to 
write down this term, and the free-field expression (B5), one finds that the overlap in (C2| is 

1 



V2 



dk_ 



V2 



dk_ 



di,e»('-+Vi)« j dze-*(«--''/^'''£:,_(0,x+)£:;_(0,a;+)e"*'=-(''-") 

-L 



(C3) 



l_ 



k_ + i/L q_ — k_ ~ i/L 



carrying out the v and z integrations to arrive at the second line. We now bring in the momentum integrals and 
change variables l_ a — L{1_ — k_) and q_ ^ b = L{q_ — k_): 



dk_ 



-ip{x2 —^i ) 



/( 



a _^ ^a/L+k_ (0, X+) ^-i(a+i)(l+x- /L) 



J^_^^JiL+kJ^:^^^{b-^)^l+x- /L) 



(C4) 
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Our task now is to take the L — > oo limit. The behaviour of the £ functions in this hmit depends crucially on the 
relative values of x'^ and fc_ , as these determine whether or not we hit the singularity. Explicitly, we have 



fa/L+fc_ (0,a;+) = exp 



dsa;^(s) 



(C5) 



We have k_ > 0. It is therefore clear that if k_ — eA_{x^) > there are no problems taking the i — >■ oo limit. 
Although we cannot carry out the s-integral exactly we know that £ becomes a phase, and this will cancel between 
the two large bracketed terms in (C4). We can therefore write part of the solution immediately: 



:0\Mx\x^)i;l{x^,x-,)\0] 



V2 



dk_ 



\f{k-)?e 



2'K a - 



27r 5- 



(C6) 



dk_ 



The integral over the range k_ € . . . e^_(a;+) requires more care. In this case there is always a value of s such that 
the integral in (C5) acquires an imaginary part which does not drop out of (|C4[). To identify the imaginary part we 



(C7) 



change variables eA_{s) = t, i.e. s = Xt and then expand the resulting t dependence in the numerator: 

eA-{x^) _ _ eA_(a;+) 



fa/L+fc_(0,2:+) = exp 



t + 



a-\-i 



exp 



iX{k_) 



dt 



where we have shown only the first term in the series and defined the function 



Hp) 



io^Xp) 



2\e\E«{Xp) ■ 

Carrying out the integrals and then taking the limit L ^ oo one finds 



£a/L+k^ (0,a;+) exp 



iX{k_){ - log[|l - eA_{x+)/k^W - in + real) 



(C8) 



(C9) 



where the real logarithm and in come from the first term of the Taylor expansion, and all other terms give real, 
a-independent contributions. (This behaviour is different from that in the second commutator term calculated in 21j, 
where the singularity always lies precisely at one of the integral's limits, rather than between them, and the a and b 
integrals are more complex.) Thus, each £ function contributes the exponential of — ttA and the rather simple, final 
expression for the first term of the density ( C2 ) is 



1 



(0|V+(a:+,X2)V^|(x+,xr)|0) = ^A 



dk_ 



\fik-re 



2 -ik-{x2 -x^ ) 



+ 



dfc_ 



i/(A-)re 



2 —ik-{x2 —x-^ ) —27rX{k-) 



eA^{x+) 

(CIO) 

Before moving on to the DHW function itself, it is worth checking this result. The second commutator term has been 
calculated independently in [2T|, see equations (4.2) and (4.9) therein. That calculation is in two dimensions, so if we 
turn off our transverse dependences completely, one finds, adding their result to our own, (CIO I, then we correctly 
obtain the (1+1) field anticommutator (Bl). Given this positive result, we return to four dimensions, and calculate 
the second term of the density U from (Bl ). The final expression for the density is, removing the test functions, 



U ■ 



-ik±y 



-A^ 



dk_ 



eA_{x+) 



Sign(fc„ — eA_{x'*'))e 



dk_ 



where we have written X2 



X + y/2, X2 = X — y/2 as in Sect, [nj and A in this expression is 

+ [k^-eA^{Xk_)Y 



A 



2\e\E\\{Xk_) 



(Cll) 



(C12) 
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We now perform the Fourier transform in ([2|. Including the Wilson line, the transformation sets 

+ eA^{x*) — , and p_ + eA_{x^) — k_ . 
The final expression for the DHW function is, over the duration of the pulse, 

W{x+;x,p) = Sign(p_) + 2P6'(-p_)0(e^_(a;+)+p_) , 
where the full expression for P is 

vrm^ + tt[p^+ eA^{x+) - eA^(Xp_+e^_(^+)))] ' 



exp 



|e|i;ii(Xp_+,^_(,+)) 



(C13) 
(CM) 

(C15) 



This is the expression given in equation (27): the more compact and intuitive notation used therein is explained in 
the text. Note that P is gauge invariant by definition, and hence for example, is actually the lightfront integral 
of the electric field, see (11). 
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